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Abstract. - In perturbative calculations of quantum-statistical zero-temperature path inte- 
grals in curvilinear coordinates one encounters Feynman diagrams involving multiple temporal 
integrals over products of distributions, which are mathematically undefined. In addition, there 
are terms proportional to powers of Dirac 5-functions at the origin coming from the measure of 
path integration. We give simple rules for integrating products of distributions in such a way 
that the results ensure coordinate independence of the path integrals. The rules are derived by 
using equations of motion and partial integration, while keeping track of certain minimal fea- 
tures originating in the unique definition of all singular integrals in 1 — e dimensions. Our rules 
yield the same results as the much more cumbersome calculations in 1 - e dimensions where 
the limit e — > is taken at the end. They also agree with the rules found in an independent 
treatment on a finite time interval. 



Introduction. - While quantum mechanical path integrals in curvilinear coordinates have 
been defined uniquely and independently of the choice of coordinates within the time-sliced 
formalism , a perturbative definition on a continuous time axis poses severe problems which 
have been solved only recently |^,^{. To exhibit the origin of the difficulties, consider the 
associated partition function calculated for periodic paths on the imaginary-time axis r: 



where A[q] is the cuclidcan action with the general form 



A[q] = / 
Jo 



dT 



■9^( q (rW(rW(r)+V( q (r)) 



(1) 



(2) 



The dots denote r-derivatives, g^iq) is a metric, and g — detg its determinant. The path 
integral is defined perturbatively as follows: The metric g^{q) and the potential V(q) are 
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expanded around some point q$ in powers of 5q^ = q^ — . After this, the action A[q] 
is separated into a free part Aa[qo;Sq] = J Q dr [^g flt/ (qo)dtSq IM dt5q v + ^^Sq^Sq"] , and an 
interacting part A in t[q ; 5q] = A[q] - Abo; Sq]. 

A first problem is encountered in the measure of functional integration in Taking 
y/ g(q) into the exponent and expanding in powers of 5q, we define an effective action A^g = 

— ^6(0) Jq dr log[g(qo + 6q)/g(qo)] which contains the infinite quantity 6(0), the (5-function at 
the origin. It is a formal representation of the inverse infinitesimal lattice spacing on the time 
axis, and is equal to the linearly divergent momentum integral J dp/(2Tt). 

The second problem arises in the expansion of Z in powers of the interaction. Performing 
all Wick contractions, Z is expressed as a sum of loop diagrams. There are interaction terms 
involving Sq 2 Sq n which lead to Feynman integrals over products of distributions. The diagrams 
contain three types of lines representing the correlation functions 

A(t-t') = (5q(r)Sq(r')} = , (3) 

d T A(r-r') = {Sq(T)6q(r')) = , (4) 

a r 9 T ,A(r-r') = (i«(r)f = (5) 

The right-hand sides define the line symbols to be used in Feynman diagrams to follow below. 
Explictly, the first correlation function reads 

A(r,r') = ^e-I— 'I. (6) 

The second correlation function has a discontinuity 

d T A(T,r') = ~e(T-r')e-^- T \ (7) 

where 

e(r - r') = -1 + 2 f dr" 5(t" - t') (8) 



is a distribution which vanishes at the origin and is equal to ±1 for positive and negative 
arguments, respectively. The third correlation function contains a <5-function: 

d T d r , A(r, r') = S(t - r') - | e -l— 'I, (9) 

Mathematically, the temporal integrals over products of such distributions are undefined jij . 
In this paper we specify these integrals by imposing the natural requirement of coordinate 
independence of the path integral (Q). By the perturbative calculation up to three loops we 
show that this requirement alone can not fix uniquely values of all ambiguous integrals over 
products of distributions. However we define a simple consistent procedure for calculating 
singular Feynman integrals. All results obtained in this way ensure coordinate independence. 
They agree with what we have obtained in our previous work. In Ref. §, we have shown 
that Feynman integrals in momentum space can be uniquely defined as e — > -limits of 
1 — e -dimensional integrals via an analytic continuation a la 't Hooft and Veltman Q . This 
definition makes path integrals coordinate- independent. In Ref. || we have given rules for 
calculating the same results directly from Feynman integrals in a 1 — e -dimensional space. 

The calculation procedure developed in this paper avoids the cumbersome evaluation of 
Feynman integrals in 1 — e dimensions. In fact, it does not require specifying any regularization 
scheme. As a fundamental byproduct, it lays the foundation for a new extension of the theory 
of distributions, in which also integrals over products arc defined, not only linear combinations. 



H. Kleinert and A. Chervyakov: Integrals over Products of Distributions. 



3 



Perturbation Expansion. The relations between singular Feynman integrals will be 
derived form the requirement of coordinate independence of the exactly solvable path integral 
of a point particle of unit mass in a harmonic potential to 2 x 2 /2, whose action is 



A, = \ I dr[x 2 (T)+u; 2 x 2 (T)] 



(10) 



For a large imaginary-time interval /?, the partition function is given by the path integral 

Z, = l'vx(r)e- A »W = e -(i/2)Trio g (-a 2 +^) = e -/^/2_ (n) 



For simplicity here and in the following the target space is assumed to be one-dimensional. 
In a higher dimensional euclidean space the problem of coordinate independence of the path 
integral (11) was first examined in the papers [|| and [Q where different noncovariant quan- 
tum corrections to the classical action ( |Io| ) were found. When this space being transformed 
holonomically to the curvilinear coordinates all covariant structures, such as the scalar curva- 
ture, remain, however, zero. Therefore the same problem can be reexamined directly in one 
dimension without loss of generality. 

A coordinate transformation turns ( Jl~T| ) into a path integral of the type (|l|) with a singular 
perturbation expansion. From our work in Refs. we know that all terms in this expan- 

sion vanish in dimensional regularization, thus ensuring the coordinate independence of the 
perturbatively defined path integral. In this paper, we proceed in the opposite direction: we 
require the vanishing of all expansion terms to find the relations for integrals over products 
of distributions. 

For simplicity we assume the coordinate transformation to preserve the symmetry x <-> —x 
of the initial oscillator, such that its power series expansion starts out like x(t) — f(q(r)) = 
q — gq 3 /3 + g 2 aq 5 /5 — • • • , where g is a smallness parameter, and a some extra parameter. 
We shall see that the identities are independent of a, such that a will merely serve to check 
the calculations. The transformation changes the partition function (Jill) into 



Vq(r) 



-Aj[q] p -A[q] 



(12) 



where is A[q] is the transformed action, whereas ^4j[g] an effective action coming from the 
Jacobian of the coordinate transformation: 



Aj[q] = -S(0) / drlog 



5q(r) ' 

The transformed action is decomposed into a free part 

Mq] = \ J dr[q 2 (r)+c 2 q 2 (r)}, 
and an interacting part, which reads to second order in g: 

Am[q] = drl -g 



(13) 



(14) 



9/ , 2 

2q 2 (r)q 2 (r) + ^q\r) 



(15) 
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To the same order in g, the Jacobian action (fi~3) is 



Aj[q} = -5(0) / dr 



-gq 2 {r) + g 2 



q 4 (r) 



(16) 



For g = 0, the transformed partition function ( |T^ ) coincides, of course, with ([TT|). When 
expanding Z of Eq. ( [l2| ) in powers of <?, we obtain Feynman integrals to each order in g, 
whose sum must vanish to ensure coordinate independence. By considering only connected 
Feynman diagrams, we study directly the ground state energy. 

Ground State Energy. - Here the coordinate independence will be tested perturbatively 
up to three loops. The graphical expansion for the ground state energy has the following 
structure: To given order g n , there exist Feynman diagrams with L = n + l,n, and n — 1 
number of loops coming from the interaction terms (jl5|) and ( fllf ) , respectively. The diagrams 
are composed of the three line types in (|^)-@, and new interaction vertices arising for each 
power of g. The diagrams coming from the Jacobian action ( [l6| ) are easily recognized by 
accompanying factors <5 n (0). 

To first order in g, there exists only three diagrams, two originated from the interaction 
(|l5|), and one from the Jacobian action ([l6|): 



9 GO 



gu: 



CO +95(0) O 



(17) 



To order g 2 , we distinguish several contributions. First there are two three-loop local diagrams 
coming from the interaction (|i~5|), and one two-loop local diagram from the Jacobian action 



(18) 



We call a diagram local if it involves no temporal integral. The Jacobian action ( Jilf ) contributes 
further the nonlocal diagrams: 



9_ 
' 2\ 



\ 1 2s 2 (o) o - 4<s(o) [ QO + CO + 2 ^ OO ] } • 



(19) 



The remaining diagrams come from the interaction ( Jl 5[ ) only. They are either of the three- 
bubble type, or of the watermelon type, each with all possible combinations of the three line 
types (|^)-(^): The sum of all three-bubbles diagrams is 



2! 



4 OX) ^ 

coo 



8w 



2 000 

s^aoo 



2 ooo 

s- 4 000] 



The watermelon-like diagrams contribute 



-^4 
2! 



(20) 



(21) 



Since the equal-time expectation value (q(r) q(r)) vanishes according to Eq. (Q) there are, in 
addition, a number of trivially vanishing diagrams, which have been omitted. 

In our previous work ^) , all integrals were calculated in D = 1 — e dimensions, taking 
the limit e — > at the end. In this way we confirmed that the sums of all Feynman diagrams 
contributing to each order in g vanish. Here we proceed in the opposite direction and derive 
the rules for integrating products of distributions from the vanishing of the sums. 
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Imposing Coordinate Independence. - In a first step we simplify the above perturbation 
expansion of the ground state energy using the inhomogeneous field equation satisfied by the 
correlation function (@): 



A(r) = - dk 



k 2 +L0 2 

and the following well-defined integrals for two correlation functions: 



(It 



A 2 (t)+u 2 A 2 {t) =A(0), 



as well as the integrals containing four correlation functions: 

jdrA\r) = ^A 3 (0), 



drA 2 (r)A 2 (r) = ±A 3 (0), 



drA 4 (r) = ^lu 2 A 3 (0). 



(22) 

(23) 

(24) 
(25) 
(26) 



These integrals are obtained directly by substituting the explicit representations (^) and ([?]) 
into (^3|), (|24|), (^5|) and (p6|), respectively. Also we shall use the singular integral obtained 
by the same direct substitution of two correlation functions: 



dr 



A"(t)+2w 2 A 2 (t)+w 4 A z (t) = / drS\T), 



4 a2/ 



(27) 



where the last integral containing the square of the <5-function will be specified later by the 
requirement of coordinate independence. 

Consider now the perturbation expansion of the ground state energy. To first order in g, 
the sum of Feynman diagrams ( |l7j ) must vanish: 



or.) + ^ 2 oo - m O=o. 

The analytic form of this relation is 

-A(0) +w 2 A(0) - 5(0)1 A (°) = °> 



(28) 



(29) 



the zero on the right-hand side being a direct consequence of the equation of motion (22) for 
the correlation function at origin. 

To order g 2 , the same equation reduces the sum of all local diagrams in ( |l8| ) to a finite 
result plus a term proportional to (5(0): 



3 (l + a) A(0) + 15 f i 



IS 



+ if u; 2 A(0) -3a 



1 



5(0) 



A 2 (0): 



35(0) --^ 2 A(0) 



Representing right-hand side diagrammatically, we obtain the identity 



E (18) = 35(0) OO " 



A 2 (0). 



(30) 



where S (frq) denotes the sum of all diagrams in Eq. (|18|). Using (t23j) together with the field 
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equation (E2J) , we reduce the sum @ of all one and two-loop bubbles diagrams to terms 
involving 6(0) and <5 2 (0): 

-|j j<5 2 (0) Jdr A 2 (t) - 2(5(0) Jdr [a(0)A 2 (t)- A(0)A 2 (t) + 2w 2 A(0)A 2 (t) 

= 25(0)A 2 (0) + S 2 (0)JdTA 2 (T). (31) 

Hence we find the diagrammatic identity 

E(19) = 2£(0) 00 + <5 2 (0)O- (32) 

Now, the terms accompanying S 2 (0) turn out to cancel similar terms coming from the sum of 
all three- loop bubbles diagrams in (|20|). In fact, the relations (^3|) and (^) lead to 



2! 



(It 



-2A(0)A(0)A 2 (t) + A 2 (0)A 2 (r) + A 2 (0)A 2 (r) + 4j 2 A 2 (0)A 2 (t) 



-4w 2 A(0)A(0)A 2 (t)+4w 4 A 2 (0)A 2 (t) = - 



dT8 2 {r) +25(0) 



A 2 (0)-5 2 (0) / drA 2 (r) 



Thus we find the diagrammatic identity for all bubbles diagrams 

S (19) + E © = - / dr5 2 (T) OO ■ 

The first two watermelon diagrams in Eq. ( pl| ) correspond to the integrals 

h = j drA 2 (r)A 2 (r), 
h = I drA(r)A 2 (r)A(r), 



(33) 

(34) 
(35) 



whose evaluation is subtle. Consider first the integral ( p4| ) which contains a square of a 
5-function. We separate this out by writing 



h = I dr A 2 (r)A 2 (r) = lf v + jf , 



with a divergent and a regular part 



/f v = A 2 (0) / drS 2 {r), 



dr A 2 (t) A 2 (t) -(5 2 (r) 



(36) 



(37) 



All other watermelon diagrams (|2lj) lead to the well-defined integrals (|26|), fl25|), and 
respectively. Substituting these and (^), (|3^) into ( ^l| ) yields the sum of all watermelon 
diagrams 



i 1 >> 



A 2 (r)A 2 (r) + 4A(t) A 2 (t)A(t) + A 4 (r) + 4w 2 A 2 (t)A 2 (t) + -w 4 A 4 (r) 

o 



= - 2A 2 (0) J dT6 2 (T)-2(l? + 47 2 ) -yw 2 A 3 (0). 



(38) 
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Adding these to (BCf), (f33J), we obtain the sum of all second-order connected diagrams 



E (all) = 3 



5(0) - / (1t8 2 (t) 



A 2 (0)-2(I? + 4I 2 )- 7 -u; 2 A 3 (0), 



(39) 



where the integrals if" and I 2 are undefined so far. This sum has to vanish to guarantee 
coordinate independence. We therfore equate to zero both the singular and finite contributions 
in Eq. (f39|). The first yields the rule for the product of two 5- functions 



S 2 ( T ) = 5(0) 6 (t) 



(40) 



This equality should of course be understood in the distributional sense, it holds after mul- 
tiplying it with an arbitrary test function and integrating over r. As it was shown in our 
previous paper the rule ( ^0|) guarantees the cancellation of all short-distance singularities 
to any order of the perturbation theory. This happens independently of the value of 5(0), 
making a regularization superfluous. There is perfect cancellation of all powers of $(0) arising 
from the expansion of the Jacobian action, which is the fundamental reason why the heuristic 
Veltman rule of setting 5(0) = is applicable everywhere without problems. 

The vanishing of the regular parts of ([39]) requires the integrals (35) and (pf) to satisfy 



I? 



4/ 2 = 



4 v ; 



7 

32^ 



(41) 



At this point we run into two difficulties. First, this single equation (|4l]) for the two un- 
defined integrals 1^ and 1% is insufficient to specify both integrals, so that the requirement 
of reparametrization invariance alone is not enough to fix all ambiguous temporal integrals 
over products of distributions. Second, and more seriously, Eq. (41) leads to conflicts with 
standard integration rules based on the use of partial integration and equation of motion. Let 
us apply these rules to the calculation of the integrals I-f and I2 in different orders. Inserting 
the equation of motion ( |22j ) into the finite part of the integral ([36]) and making use of the 
regular integral (0), we find immediately 



(42) 



Jf = J cLtA 2 (t) \ A 2 (t) -5 2 (t) 

= -2^A 3 (0)+ W 4 |rfrA 4 (r) = -^ 2 A 3 (0) = -^-. 



The same substitution of the equation of motion (|2^ ) into the other ambiguous integral I2 of 
( |35l) leads, after performing the the regular integral (||]), to 



h 



dr A 2 (t) A(t)5(t) +lu 2 J dr A 2 (t) A 2 (t) 



where / denotes the undefined integral over a product of distributions 

1=1 dre 2 (T)5(T). 

This integral can apparently be fixed by applying partial integration to the integral 
reduces it to the completely regular form (Eq): 



(43) 



(44) 



5h which 



drA(r)— A 3 (r) 



drA\r) = --u 2 A\0) 



1 

96^' 



(45) 
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There are no boundary terms due to the exponential vanishing at infinity of all functions 
involved. From ( [43| ) and (45) we conclude that I = 1/3. 

At this point we observe a conflict. The results ([l5]) and ( |42| ) do not obey the equation 
( |4T| ) , such that they are incompatible with the necessary coordinate independence of the path 
integral. This was the reason to add a noncovariant quantum correction term AV = —g 2 (q 2 /6) 
to the classical action ( |To| ) in the previous paper j?j . For perturbative calculation on a finite- 
time interval it was also done in ||. 

From the perspective of our previous papers |^||] where all integrals were defined in d = 1— e 
dimensions and continued to e — > at the end, the above-observed inconsistency is obvious: 
Arbitrary application of partial integration and equation of motion to one-dimensional inte- 
grals is forbidden whenever several dots can correspond to different contractions of partial 
derivatives d a , dp, . . ., from which they arise in the limit d — ► 1. The different contractions 
may lead to different integrals. In the pure one-dimensional calculation of the integrals if 1 and 
1% this ambiguity can be accounted for by using partial integration and equation of motion 
(^2| ) only according to the following integration rules: 

1 . We perform a partial integration which allows us to apply subsequently the equation of 
motion (53). 

2. If the equation of motion (22) leads to integrals of the type (p4|), they must be performed 
using the naively the Dirac rule for the <5-function and the property e(0) = 0. 

3. The above procedure leaves in general singular integrals, which must be treated once 
more with the same rules. 

Let us show that calculating the integrals I± and li with these rules is consistent with 
the coordinate independence condition (f4l|). In the integral I2 of (|3^) we first apply partial 
integration to find 

h = ~ JdrA(r)A(r)^-[A 2 (r) 

= -I fdrA^r)- 1 - fdrA(r)A 2 (r)A(r). 



(46) 



is 



with no contributions from the boundary terms. Note that the partial integration ( |45[ 
forbidden since it does not allow for a subsequent application of the equation of motion (|22|) . 
On the right-hand side of ( filf ) it can be applied. This leads to a combination of two regular 
integrals ( p5j ) and (26) and the singular integral /, which we evaluate with the naive Dirac 
rule to / — 0, resulting in 



h = -- I drA\r) 



I - -lo z A 6 (Q) = -- 



\ J drA 2 (r)A(r)5(r)-^ 
1 



drA 2 (r) A 2 (r) 



(47) 



If we calculate the finite part /-p of the integral ( |36| ) with the new rules we obtain a result 
different from p2]). Integrating the first term in brackets by parts and using the equation of 
motion (B2fl, we obtain 



I? 



drA 2 (r) A 2 (r) - 5 2 {t) 



dr 



A (r) A(r) A 2 (r) - 2A(r) A 2 (r) A(r) - A 2 (r) 5 2 {t) 



= I dr A(t)A 2 (t)S{t)-A 2 {t)S 2 (t) - 2I 2 - u 2 I drA 2 (r) A 2 (r) . (48) 



H. Kleinert and A. Chervyakov: Integrals over Products of Distributions. 



9 



The last two terms are already known, while the remaining singular integral in brackets must 
be subjected once more to the same treatment. It is integrated by parts so that the equation 
of motion (E3) can be applied to obtain 



A(r) A 2 (t)6(t) - A 2 (t)S 2 (t) 



(It 

= -J dT 

Inserting this into Eq 



A(r) A z {t) + 2A 2 {t) A(t) S(t) - I drA^(r) 5\t) = -J* A 3 (0) - f- 1. (49) 



yields 



drA 2 (r) A 2 (t)-<5 2 (t) 



-2J a 



the right-hand side following for 1 = 0. 

We see now that the integrals p^ ) and ( |50| ) calculated with the new rules obey the equation 
( f4l|) which guarantees coordinate independence of the path integral. 

The applicability of the rules (1-3) follows immediately from the previously established 
dimensional continuation [^[^1. It avoids completely the cumbersome calculations in 1 — e- 
dimension with the subsequent limit e — > 0. Only some intermediate steps of the derivation 
require keeping track of the (/-dimensional origin of the rules. For this, we continue the 



imaginary time coordinate r to a d-dimensional spacetime vector r 



(r,x l , 



„d-l\ 



and note that the equation of motion (|22|) becomes a scalar field equation of the Klein-Gordon 
type 

{-d 2 a +lo 2 ) A(x) =tW(x). (51) 

In d dimensions, the relevant second-order diagrams are obtained by decomposing the har- 
monic expectation value 



d d x < q 2 a (x)q 2 (x)q 2 (0)q 2 (0)> 



(52) 



into a sum of products of four two-point correlation functions according to the Wick rule. The 
fields q a (x) are the ci-dimensional extensions q a {x) = d a q(x) of q(x). Now the d-dimensional 
integrals, corresponding to the integrals (|3~i| ) and (p5|), are defined uniquely by the contractions 



If 



1 I 



d d x <q a (x)qa(x)q{x)q(x)qp(Q)qp(0)q(0)q(0) > 
d d xA 2 (x)A 2 a0 (x), 



d d x <q a (x)q a (x)q(x)q(x)qp(0)qp(0)q(0)q(0) > 
d d x A(x) A a (x) Ap (x) A a!3 (x) . 



(53) 



(54) 



The different derivatives d a dp acting on A(x) prevent us from applying the field equation 
(|5l"|). This obstacle was hidden in the one-dimensional formulation. It can be overcome by a 
partial integration. Starting with I d , we obtain 



TO, _ _ 
1 9 — 



\ J d d xA 2 (x) [A 2 a (x) + A(x) A aa (x)] , 



(55) 
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Treating if likewise we find 

If = -21$ + J d d xA 2 {x)A 2 aa (x)+2 J d d xA{x)A 2 {x)A aa {x). (56) 

In the second equation we have used the fact that <9 Q A Qj a = dpA\\. The right-hand sides 
of (H|) and ( |56|) contain now the contracted derivatives d 2 such that we can apply the field 
equation fell). This mechanism works to all orders in the perturbation expansion which is 



the reason for the applicability of the rules (1) and (2) which led to the results (47) and (|30j) 
ensuring coordinate independence. 

The rule (3) is a consequence of regularized equation in d = 1 — e dimension 

A(x) A 2 a (x) 5(x) ~ \x\ 2e 5{x) = , (57) 

which makes the singular product of distributions in Eq. ( ff4| ) vanishes before taking e — > 0- 
limit. This rule is in agreement with the earlier result derived for finite time intervals in 
Rcf. g. 

Let us illustrate how different contractions of partial derivatives d a , dp, . . ., may lead to 
different integrals in the d — ■+ 1-limit. The simplest example is the anomalous integral 

It = J d d xA 2 {x)A 2 a(3 {x) - J d d xA 2 (x)A 2 aa (x), (58) 

where the integrals on the right-hand side arc indistinguishable in d — 1 dimension. However, 
it follows immediately from Eqs. (|55|) and (|5^ ) that 



I d n = ld d xA 2 a (x)A 2 (x) + 3jd d xA(x)A 2 (x)A aa (x) 

d d x A 2 a (x) A} (x) + 'Slu 2 J d d x A 2 (x) A| (x) 

- 3 J d d x A(x) A 2 (x)6(x) d -+ i-. (59) 

The dimensional continuation of all other second-order diagrams does not change the results 
of one-dimensional calculation. Only the second diagram in Eq. ( pcf ) requires taking more 
care. However, the corresponding <i-dimensional integral has the property 

A 2 (0) / d d xA 2 (x) = A 2 (0) / d d xA 2 aa (x) , (60) 



which allows us to apply the equation of motion ( pl|) directly as in one dimension. Thus, 
by keeping only track of a few essential properties of the theory in d dimensions we indeed 
obtain a simple consistent procedure for calculating singular Feynman integrals. All results 
obtained in this way ensure coordinate independence. Our procedure gives us unique rules 
telling us where we are allowed to apply partial integration and the equation of motion in one- 
dimensional expressions. Ultimately, all integrals are brought to a regular form, which can 
be continued back to one time dimension for a direct evaluation. This procedure is obviously 
much simpler than the previous explicit calculations in <i-dimension with the limit d — * 1 
taken at the end. 

The above rules for integrating products of distributions are in complete agreement with 
the rules derived earlier in a different way for path integrals on a finite time interval where 
the infrared regulator uj 2 q 2 was not needed in the free part of the action Aq- 
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Let us briefly discuss the alternative possibility of giving up partial integration completely 
in ambiguous integrals containing e- and 5-function, or their time derivatives, which makes 
unnecessary to satisfy Eq. ( [45| ) . This yields a freedom in the definition of integral over product 
of distribution which can be used to fix I = 1/4 from the requirement of coordinate 
independence JlOfl . Indeed, this value of / makes the integral ( fl3| ) equal to I2 = such 
that (^) is satisfied and coordinate independence ensured. In contrast, giving up partial 
integration, the authors of Refs. |3 
the integral / should vanish as well: 1 = 0. Then Eq. ( [43| ) yields I2 = l/32w which together 
with ( |42| ) does not obey the coordinate independence condition (f5|), making yet an another 
noncovariant quantum correction AF = g 2 (q 2 /2) necessary in the action, which we reject 
since it contradicts Feynman's original rules of path integration. We do not consider giving 
up partial integration as an attractive option since it is an important tool for calculating 
higher- loop diagrams. 

Summary. - We have defined singular Feynman integrals in perturbative calculation of 
zero-temperature path integrals in such a way that coordinate independence is guaranteed. To 
second-order in perturbation theory, one encounters two singular integrals containing products 
of distributions and finds one relation between them due to coordinate independence. The 
evaluation of the individual integrals was made unique by analytic regularization In this 

paper we have shown that there are simple rules for obtaining the same results without the 
analytic continuation of Feynman integrals to d dimensions. We have merely kept track of the 
different contractions of the derivatives and performed partial integrations until one obtains 
Laplace operators which allow us to use the d-dimensional field equation which removes one 
singularity. If repeated recursively, this procedure leads to regular integrals. There is no need 
to specify a regularization scheme, and calculations are much simpler than the previous ones 
in (i-dimension with the limit d — > 1 taken at the end. Our results are agreement with the 
earlier result derived for finite time intervals in Ref. Q . 

Just as in the time-sliced definition of path integrals in curved space in Ref. jlj, there is 
absolutely no need for extra compensating potential terms found necessary in the treatments 
in Refs. 1,0,0. 



, 1 1 1 have assumed the vanishing e 2 (r) at t = so that 
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